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Foreword 

This  is  the  final  report  of  AHQ  Grant  DAAG29-77-G-0051.  From  the 
beginning  of  1977  to  September  1978,  the  principal  investigator  was 
supported  20%  by  the  grant  in  the  academic  year  and  100%  in  the  summer. 

On  the  basics  of  the  work  performed  for  the  Grant  the  following 
personnel  received  degrees: 

J.  A.  Holzer,  Ph.D. 

W.  B.  Fowler,  Ph.D. 

W.  D.  Eberhardt,  M.Sc. 

Most  of  the  work  reported  here  was  performed  in  1979  and  1978. 

The  principal  investigator  was  on  his  sabbatical  from  October  1978  to 
June  1979.  During  this  period  he  was  not  supported  by  the  Grant  and 
no  effort  was  made  for  the  task  except  writing  up  the  results  for 
publication.  The  report  is  written  in  June  1979,  when  he  was  away 
from  the  Oak  Ridge  National  Laboratory  where  he  held  a summer  appoint- 
ment. 

In  Section  1 the  accomplishment  of  the  works  performed  under  the 
support  of  the  grant  is  summarized  in  three  categories.  The  first  two 
categories  are  directly  related  to  the  original  proposal  and  are 
discussed  in  details  in  Section  2A  and  2B.  The  last  part  contains  a 
summary  of  the  research  in  some  other  area.  The  abstract  of  all 
the  publications  is  given  in  Appendix  A,  An  important  part  of  the 
effort  in  the  speckle  statistics  deals  with  the  computation  of  the 
probability  density  function  of  light  received  by  a finite  size 
detector.  This  work  is  included  as  Appendix  B to  be  submitted  for 
publication  later. 
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Section  1 - Outline  of  the  work  performed 

The  following  is  a list  of  publications  of  the  principal  investi- 
tor  arranged  in  chronical  order  after  the  grant  was  started  in  the 
beginning  of  1977. 

A 1.  "Scattering  of  Electromagnetic  Waves  from  a Rough  Surface  II" 
with  Holzer,  J.  of  Appl.  Fhys.  49,  1002  (1978). 

A 2.  "Scattering  of  an  Electromagnetic  Wave  from  a Very  Rough  Semi- 
Infinite  Dielectric  Plane  (Exact  Treatment  of  the  Boundary 
Condition)"  J.  of  Appl.  Phys.  49,  994  (1978)  with  W.  D.  Eberhardt. 

3.  "Magnetic  Phase  Transitions  in  the  Dlcke  Model,"  C.  Bowden,  J. 
of  Phys.  A 11,  151  (1978). 

4.  "Distorted  Wave  Born  Approximation  for  the  Rate  Constant  of  V-V 
Energy  Transfer,"  with  Stettler  and  Wltriol,  J.  of  Chem.  Phys. 

69,  3112  (1978). 

A 5.  "Explanation  of  the  Experimental  Results  of  Light  Back  Scattered 
from  a Rough  Surface,"  with  W.  D.  Eberhardt,  J.  of  Opt.  Soc. 

(AM)  68,  323  (1978). 

C 6.  "Radiative  Energy  Transfer  through  Fog"  Appl.  Optics  17,  1997  (1978). 

A 7.  "Scattering  of  Light  from  Rayleigh  Waves"  Appl.  Phys.  Lett.  32, 

399  (1978). 

C 8.  "Cooperative  Behavior  among  Three  Level  System  I:  Transient  Effects 
of  Coherent  Optical  Pumping,"  Phys.  Rev.  18  A 1558  (1978)  with 
C.  Bowden. 
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B 9.  "Probability  Density  Function  of  Scattered  Light  from  a Moderately 
Rough  Surface,"  Appl.  Opt.  18,  312  (1979  with  Stettler  and  Holzer. 
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C 10.  "Scattering  of  Light  from  a Dielectric  Body  of  Irregular  Shapes," 
J.  of  Opt.  Soc.  69,  756  (1979)  with  B.  Fowler. 

C 11.  "Phase  Transition  In  the  Three-Level  Dicke  Model",  J.  of  Phys.  A 
(1979)  with  C.  Bowden. 

B 12.  "Statistical  Properties  of  Reflected  Light  from  a Moderately 

Rough  Surface  Observed  through  Finite-Size  Aperature",  this  work 
is  included  in  the  final  report  as  Appendix  B. 


All  the  papers  are  published  in  refereed  journals.  Papers 
delivered  in  meetings  and  conferences  are  not  included. 

Publication  3 and  4 are  the  results  of  much  earlier  work  and  will  not 
be  discussed.  We  divide  all  other  publications  in  three  categories. 
Category  A and  B deal  with  light  reflection  from  a rough  surface  and 
speckle  statistics  that  are  the  main  topics  which  were  proposed  to  be 
studied  in  the  proposal. 


) 


In  publication  A we  show  the  limitation  of  the  conventional  ap- 
proach of  the  scattering  of  light  from  a rough  surface  by  using  the 
vectorial  Kirchhoff  integral  based  on  Beckmann  and  Splzzichino. ^ It 
is  clear  that  this  approach  is  good  for  a good  conductor  and  poor  for 
a dielectric  body.  A realistic  surface  model  is  constructed  to  explain 
some  experimental  data  and  the  importance  of  the  boundary  condition 
is  demonstrated. 

Since  the  speckle  statistics  of  light  reflected  from  a rough 
surface  has  been  discussed  exclusively  for  a very  rough  surface,  we 
compute  the  probability  density  function  P(I)  for  a moderately  smooth 
surface  both  for  a point  end  finite  site  detector.  We  show  how  P(I) 
can  be  expressed  as  a function  of  the  surface  parameters. 

During  1977  and  1978,  we  also  worked  on  the  laser  propagation  in 
various  mediums  as  a part  of  the  Ph.D.  dissertation  of  Mr.  Fowler. 

There  was  also  active  collaboration  with  Dr.  C.  Bowden  on  quantum 
electronics.  A brief  discussion  of  these  works  and  the  various 
contributions  are  given  in  the  section. 
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Section  2 

A.  Light  Scattering  from  a Rough  Surface 

Most  works  in  the  analysis  of  experimental  data  of  light 

(0 

scattering  of  a very  rough  surface  use  Beckmann  and  Spizzlchino 's  work 

2-4 

as  standard  reference.  In  much  earlier  works,  we  have  shown  that 

the  so-called  satisfactory  results  based  on  Ref.  1 is  actually  a 

combination  of  many  approximations,  both  numerical  and  analytical 

involving  the  treatment  of  random  variables  . 

In  Ref.  5,  this  issue  was  carefully  studied.  It  is  shown  that  by 

using  the  approximate  boundary  conditions  to  relate  the  scattered 

field  amplitude  E and  the  incident  field  amplitude  E. 

sc  in 

E,c  ■ (1  ♦ r>  hn  m 

where  the  sign  is  chosen  for  different  polarization  is  a very  poor 
approximation  for  a dielectric  body.  This  is  illustrated  in  the 
following  figure. 
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Here  -the  results  for  very  large  dielectric  constant  $ are  almost 
Independent  of  the  boundary  conditions,  whereas  for  0,  Eq.  (1) 
predicts  a reflection  of  finite  magnitude  although  It  Is  clear  that 
no  light  should  be  reflected  from  a non-existent  boundary.  It  should 
be  emphasized  that  this  conclusion  Is  reached  without  using  any 
questionable  approximation  In  the  treatment  of  random  variables  and 
numerical  integration. 

The  vectorial  Kirchhoff  integral  used  in  Ref.  1 also  fails  to 

predict  depolarization  of  back-scattered  light  from  a rough  surface, 

contradicting  the  well-known  experimental  results.**  In  Ref.  7 and  8, 

a realistic  surface  model  Is  constructed,  so  that  the  surface  becomes 

a combination  of  independent  cells  which  are  described  by  a quadratic 

equation.  In  other  words,  this  surface  generalizes  the  tangent  plane 

approximation  to  include  the  curvature  correction.  When  a calculation 

7 8 

including  exact  treatment  of  the  boundary  conditions  is  completed,  ' 
we  find  that  the  results  not  only  predict  depolarization  but  also  have 
correct  angular  distribution.  In  the  following  figure,  the  dashed  line 
is  roughly  the  experimental vresult,  and  the  solid  lines  are  from  the 
calculation.  Notice  that  the  magnitude  of  the  depolarized  component  is 
not  sensitive  to  the  surface  roughness  in  contrast  to  the  co-polarlzed 
component.  This  result  is  also  in  good  agreement  with  the  observation. ^ 
We  are  not  aware  of  any  work  in  the  literature  based  on  the  first 
principle  calculation  that  can  derive  these  results. 


•or  dMforant  InddM  anglw  t and  4.  ft  “ 0.1.  < - 2.6. 
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The  only  undesirable  feature  Is  that  only  small  depolarization 
is  predicted  whereas  the  experimental  value  is  much  larger.  It  is 
possible  that  this  deficiency  is  cause?,  by  neglecting  multiple 
scattering  in  the  calculation. 

The  importance  of  the  boundary  conditions  in  any  calculation  of 
the  scattering  of  light  is  again  demonstrated  in  Ref.  9 , where  the 
laser  light  reflection  from  a corrugated  surface  generated  by  acoustic 
waves  propagating  in  a dielectric  medium  is  computed.  Experimental 
data  and  previous  works  on  this  subject  are  reviewed  by  Lean.1® 

We  show  that  previous  calculation  based  on  Ref.l  actually  never  predict 
even  the  shift  of  the  Brewster  angle  for  diffracted  wave  of  different 
orders.  This  shift  has  been  measured  experimentally  but  unfortunate 
algebraic  errors  in  the  earlier  work1®  give  a wrong  impression,  that  the 
calculation  based  on  Ref.  1 is  satisfactory.  Our  calculation  satisfies 
the  boundary  conditions  exactly  and  should  be  useful  for  future 
experimental  work  in  this  area. 

11-12 

We  should  mention  many  works  in  the  recent  literature  on  the 

light  scattering  from  a moderately  rough  surface,  when  the  perturbation 
technique  is  used  to  satisfy  the  boundary  conditions.  Our  works,  as  a 
matter  of  fact,  follow  their  technique  closely,  with  the  important 
difference  that  we  deal  with  very  rough  surfaces  where  considerable 
modifications  are  required  as  we  did  throughout  these  publications. 
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II  B.  Speckle  Statistics 

The  reflected  intensity  from  a rough  surface  as  discussed  in 

II  A is  related  to  the  surface  parameters.  More  information  on  the 

characteristics  of  the  surface  can  be  obtained  by  measuring  the 

probability  density  function  P(I),  which  is  determined  by  scanning 

across  the  illuminated  rough  surface.  This  subject  has  been  reviewed 

13  14 

In  the  literature  by  Goodman  and  Dainty  who  provided  a complete 
description  of  the  background  material.  There  are  also  a great  number 
of  experimental  activities. ^ ^ 

In  both  the  earlier  works  and  later  ones,  P(I)  is  computed  by 
assuming  Gaussian  statistical  properties  of  the  scattered  amplitude; 
real  and  imaginary  part  of  the  amplitude  from  a circular  Gaussian 
variable.  This  approach  is  similar  to  those  in  the  random  signed 
analysis  The  coherent  component  is  considered  as  the  signal  and  the 
incoherent  component  the  noise,  and  is  very  successful  for  light  reflec- 
ted from  a very  rough  surface,  since  the  phase  scattering  amplitude 
should  vary  between  zero  and  2tt  uniformly.  The  calculation  of  P(I)  for 
this  case  seems  to  be  complete. 

In  the  past  year,  we  concentrated  on  P(I)  for  reflected  light  from 
a moderately  rough  surface  a £ X when  a is  the  rms  of  the  roughness 
and  X is  the  wavelength.  In  this  case  the  problem  is  somewhat 
simplified  since  the  amplitude  can  be  expanded  in  o/X  . On  the  other 
hand,  the  statistical  properties  about  the  amplitude  assumed  in  the 
previous  works  are  not  valid  and  P(I)  should  be  expressed  in  terms  of 
the  parameters  of  the  surface  characteristics.  This  is  what  we  have 
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accomplished  In  Refs. 16  and  17.  In  Che  first  paper:  P(I)  for  a point 
detector  Is  considered  l.e.  the  size  of  the  detector  Is  less  than  that 
of  one  speckle.  P(I)  Is  derived  without  using  any  specific  assumption 
about  the  surface  properties  except  the  Illuminated  area  must  be  much 
greater  than  the  size  of  the  correlation  cell  defined  In  the  surface 
model.  This  assumption  is  sufficient  for  us  to  invoke  the  Central 
Limit  theorem  to  obtain  P(I)  In  terms  of  the  first  two  moments  of  the 
surface  variables  Z(x).  The  final  expression  for  P(I)  contains  three 
Independent  parameters  related  to  the  rough  surface.  We  have  used 
this  expression  to  fit  the  data  by  U.  S.  Missile  Command,  Huntsville, 
Alabama  with  Improved  accuracy. 

In  Ref.  17,  which  is  included  as  Appendix  B,  we  generalize  the 
work  in  Ref. 16  for  detector  of  finite  size,  where  several  speckles 
may  appear  and  one  has  to  consider  the  fluctuation  of  many  speckles. 
Again  this  problem  differs  from  our  case  to  the  very  rough  surface  case 
in  that  we  have  only  one  real  random  variable  Z(x)  instead  of  the  real 
and  imaginary  component  of  the  amplitude.  For  a complete  general  case 
of  arbitrary  scattering  geometry,  the  numerical  work  is  too  formidable, 
as  indicated  in  Ref.  17.  We  work  out  for  a case  of  practical  interest 
and  the  numerical  results  are  expressed  in  terms  of  relevant  parameters, 
average  number  N of  speckles  in  the  detector,  (o/X)  , and  the  ratio 
of  correlation  length  to  the  size  of  the  illuminated  area.  It  is  shown 
that  for  N » 1,  P(I)  is  not  remarkably  different  from  that  of  a very 
rough  surface,  if  the  parameters  are  properly  identified.  For  N ^ 1 
a careful  numerical  solution  is  necessary  for  accurate  results. 
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II  C.  Miscellaneous 

RefsJ.8-21  Che  product  of  other  research  activities  beside  the 

grant,  although  they  are  all  related  in  a broad  sense.  In  Ref.  18, 

the  radiative  transfer  equation  is  solved  numerically  for  a medium 

consisting  of  aerosols  for  propagation  of  a coherent  light  source  of 

arbitrary  Intensity  distribution.  The  contribution  is  mainly  numerical, 

but  it  should  be  useful  for  practical  simulation,  and  important  to 

22 

the  application  of  laser  radar  technology,  since  no  limitation  is 
imposed  on  the  form  of  the  intensity  distribution  in  our  calculation. 

The  loss  of  coherence  and  contrast  as  a function  of  the  propagation 
distance  through  the  medium  is  clearly  demonstrated  in  the  numerical 
calculation. 

Continuing  the  effort  in  Ref.  18,  we  later  address  to  the  scattering 

of  light  from  a dielectic  body  of  arbitrary  size  and  shpe.  This  problem 

22 

is  of  interest  both  to  laser  radar  technology  where  the  contribution 

of  anisotropic  particles  is  a subject  of  concern  and  applied  physics 

where  the  Rayleigh's  Hypothesis  has  always  been  a subject  of  study. 

In  Ref.  19,  we  compute  the  diffracted  intensity  of  an  Imperfect  sphere 

24 

and  compare  with  the  known  exact  solution  to  show  the  limitation  of 
Rayleigh's  Hypothesis.  We  also  obtain  numerical  results  consistent 

with  the  experimental  results  in  the  light  scattering  of  an  anisotropic 

25  25 

medium.  It  is  confirmed  that  although  the  previous  empirical 

treatments  are  simple  and  efficient,  in  some  specie  cases,  their 

results,  however,  should  be  used  cautiously. 


In  the  later  part  of  1978,  the  Interests  of  the  principal 
investigator  was  shifted  toward  quantum  electronics.  Refs. 20  and  21 
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specifically  are  addressed  to  Che  three-level  molecules,  In  contrast 

26 

to  the  two-level  molecules  used  In  the  literature.  Our  model  of  course 
Is  more  realistic.  It  is  found  in  Ref.  20,  however,  that  the  character- 
istics of  the  phase  transition  in  the  Dicke  model  is  not  changed  very 
much  except  there  exist  two  second  order  phase  transitions  instead  of 
one.  One  important  result  derived  in  Ref.  20  is  that  the  transition 

3 

temperature  should  be  given  in  rms  terms  of  atoms  inside  a volume  X 
instead  of  the  V,  the  quantization  volume  which  is  arbitrary  in  the 
Dicke  model. 

In  Ref.  21,  we  also  study  the  cooperative  behavior  in  the  three- 

level  molecules.  This  model  makes  it  possible  to  compute  effect  of 

coherent  pumping.  The  results  seem  to  explain  some  of  the  well-known 

26—27 

discrepancy  between  the  two- level  model  and  experimental  data. 
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A general  aipraaainu  for  the  probability  deneity  function  P(D  of  the  intenaity  of  light  reflected  from  a mod- 
erately rough  surface  in  the  far  Raid  is  derived  by  the  cumulant  expansion  of  the  random  variable.  It  is 
shown  that  P(D  satisfies  the  modified  Rician  statistics  for  soma  general  surface  models  in  the  limit  of  large 
LIT,  where  L ia  the  dimension  of  the  illuminating  ana,  and  T ia  the  correlation  length. 
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Croeseectlopa  far  an  electromagnetic  wave  scattered  from  dislectric  bodies  of  trrwgular  ahape  have 
bean  obtained.  Urn  boundary  conditions  are  reduced  to  a sat  of  liaanr  algebraic  equations  whose 
vsriaMm  amtheaxpanaitm  cosfflcieoti  of  thaelactric^and  magnetic  finite.  Numerical  rwauhs  are 

Raylnlgh  hypothesis  implicitly  used  In  the  formulation.  The  procedures  uawd  byChyfaK  (hum  and 
Planish  to  fit  the  experimental  data  wwafaonuminad  and  diamond. 
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TT*r3-D  time-independent  unpolarixed  radiative  transfer  equation  it  developed  into  a form  amanabla  to  nu- 
maneJ  pulton  by  «P«uion  of  the  totol  intonmty  in  order,  of  acattering.  A 2-D  numerical  algorithm  i. 
ormulated  using  an  eight-direction  quadrature  approximation  to  tha  tource  function  and  an  optimized 

pat^by-prth  intensity  integrmtioo  scheme.  Computer  coda  generated  from  the  algorithm  itueedto 

ma  th,  multiple  ecrttermg  ^fccu  at  1.06  ,m.  3.0  wand  10.6*0.  in  Drinnendjian  C3  fog  on  a finite  uniform 
■ntMHty  profile  beam  and  an  infinite  Gauteian  intensity  profile  beam.  Tranamiaaion  and  backtcatter  for 
both  beams  are  examined  at  well  at  conUatt  degradation  of  the  uniform  beam  and  tpread  of  the  r- 
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Phase  Transition  In  the  Multi-Mode  Two-  and  Three-Level  Dlcke  Model 
(Green's  Function  Method) 
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Abstract.  By  using  the  thermodynamical  Green's  function  method, 
we  study  the  phase  transition  of  a collection  of  two-  as  well  as 
three-level  atoms  interacting  with  the  electromagnetic  field  and 
contained  within  a volume  much  smaller  than  the  smallest  resonance 
wavelength  (Dicke  model).  We  show  for  the  case  of  three-level 
atoms,  the  existence  of  two  critical  temperatures  where  a second 
order  phase  transition  takes  place.  The  upper  critical  tempera- 
ture Is  determined  by  the  largest  of  the  three  coupling  constants 
In  the  three-level  model,  whereas  the  lower  critical  temperature 
Is  a function  of  the  remaining  two  coupling  constants.  We  also 
show  that  for  a collection  of  two-  as  well  as  three-level  atoms, 
the  critical  temperature  depends  upon  the  number  of  atoms  contain- 
ed  within  X*  (Xr  Is  a resonance  wavelength)  rather  than  the  den- 
sity *s  has  previously  been  suggested.  Also,  using  the  Green's 
•• 

« 

function  method,  we  show  the  formal  equivalence  of  the  Dlcke  model 
• * ••» 

In 'thermodynamic  equilibrium  with  the  BCS  model  of  superconductivity. 
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The  probability  density  function  P(I)  of  the  inte- 
grated intensity  I , reflected  from  a moderately  rough 

2 

surface,  is  computed  correct  to  (2ira/X)  where  a is  the 
rms  of  the  rough  surface  z(x)  and  X is  the  incident  wave- 
length. By  using  the  assumption  that  z(x)  is  a Gaussian 
random  variable,  we  obtain  P(i)  as  a function  of  the  size 
of  the  detector  aperture,  and  the  dependence  of  the  corre- 
lation length  of  the  surface  model.  Numerical  results  and 
a comparison  with  previous  works  are  given. 
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I.  INTRODUCTION 


Statistical  properties  of  speckle  patterns  have  been  discussed 
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extensively  in  recent  literature.  Several  review  articles  provide 
excellent  background  material.  The  calculation  of  the  probability  den- 
sity function  P(I)  for  reflected  light  from  a very  rough  surface,  or 
transmitted  through  a very  diffusive  medium.  Is  based  on  the  justifi- 
able assumption  that  the  field  amplitudes  obey  Gaussian  statistics.  The 
result,  which  Is  Independent  of  the  details  of  the  surface  properties, 
is  that  P(I)  Is  a Gamma  distribution  function  with  the  order  parameter 

depending  on  the  size  of  the  detector  aperture.  Numerical  results  and 

5 6 7 

other  applications  are  given  by  Seribot  , Dainty  , and  Barakat  . 

8 9 

For  a moderately  rought  surface,  P(I)  has  been  calculated  ’ and 

experimentally^  observed  for  a point  detector.  In  References  8-10, 

P(I)  Is  calculated  by  assuming  that  the  field  amplitudes  obey  non-circular 
«.► 

Gaussian  statistics.  In  Reference  11,  we  have  transferred  the  assumption 
of  Gaussianlty  from  the  field  onto  the  surface  variable  z(x).  Since  the 
field  amplitude  and  z(x)  are  proportional  for  a moderately  rough  surface 
where  o/X  < < 1,  the  difference  between  Reference  11  and  References  8-10 
is  only  In  the  definition  of  various  moments  that  appear  In  P(I). 

We  are  not  aware  of  any  work  on  P(I)  from  a moderately  rough  surface 
for  finite-size  detector  aperture.  An  obvious  and  simple  method  is  to 
extend  the  work  by  Ohtsubo  and  Asakura  who  obtained  P(I)  for  a sum  of 
N independent  speckles.  Their  procedure  Is  similar  to  the  earlier  work 
by  Goodman2  who  derived  a Gamma  distribution  for  the  P(I)  of  a very  rough 
surface  for  a finite-size  detector.  If  we  follow  this  method,  the  derived 
P(I)  however,  contains  a parameter  N , the  number  of  speckles  In  the 
finite-size  detector.  N Is  proportional  to  the  size  of  the  aperture 
for  larger  N . For  small  N , P(I)  is  more  complicated  as  demonstrated 
by  Seribot.5  „ 


The  purpose  of  this  paper  is  to  calculate  P(I)  of  a moderately  rough 
surface  observed  through  a finite-size  aperture  correct  to  the  second  order 
in  (ct/X)  by  assuming  the  Gaussianlty  of  z(x).  Besides  this  assumption  and 
other  simplification  Involving  scattering  geometry,  no  other  assumptions 
are  made  and  we  believe  that  this  work  and  Reference  5 complete  the  cal- 
culation of  P(I)  for  surfaces  of  moderate  and  extreme  roughness. 

In  the  second  section,  we  present  the  general  principle  how  P(I)  can 
be  computed  from  the  known  statistical  properties.  In  the  third  section 
a simpler  scattering  geometry  is  considered  and  solved  numerically  in  the 
last  section. 

II.  GENERAL 

We  consider  the  field  amplitude  E(x,u)  scattered  from  a rough,  illumi- 

^ A A 

nated  surface  (x  * xi  + yj)  located  at  a distance  d from  the  detecting 

^ A A 

plane,  (u  * ui  + vj  , and  all  vectors  In  this  paper  are  two-dimensional.) 

It  Is  assumed  that  the  incident  light  is  linearly  polarized  and  the  depol- 
arization of  the  scattered  light  is  not  important  so  that  E(x,u)  is  a scalar 
function.  Generalization  of  the  calculation  to  a vectorial  E(x,u)  has  been 
discussed  by  Goodman.  We  limit  our  calculation  to  the  second  order  in 
2/X(  X*  2ir/k)  and  expand 

E(x,u)  - f0  + f,  kz  + f2k2z2  (1) 

where  f0,  fj  and  f2  are  assumed  to  be  real  functions  of  x and  u . 

The  Imaginary  part  of  E(x,u)  will  be  Included  later.  Equation  (1)  then  Is 
a general  form  for  all  scattering  geometries  from  which  f0»  fj  and  f2  are 
determined.  Our  problem  now  Is  to  derive  P(I)  for  a given  statistical 
property  of  z(x).  We  assume  that  z(x)  are  Gaussian  random  variables 
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whose  first  two  moments  are  given  by 

< z(x)  > » 0 , < z(x)  z(x* ) > * a2  exp[-  (x-x')2/t2J.  (2) 

Equation  (2)  replaces  the  assumption  on  the  Gausslanlty  of  the  E(x,u)  In  the 

2-4 

treatment  of  scattering  amplitude  of  a very  rough  surface.  Physically 
E(x,u)  Is  a macroscoplcal  quantity  whereas  z(x)  are  the  characteristics  of 
the  surface  property.  The  assumptions  used  In  this  work  and  previous  work 
are  different  In  nature.  Equation  (1),  however,  connects  these  two  quantities 
and  It  Is  clear  why  they  serve  the  same  purpose  In  the  derivation  of  P(I). 

It  should  be  emphasized  that  because  of  Eq.  (2),  P(I)  In  this  work  depends 
only  on  the  surface  properties. 

The  major  difference  between  our  problem  defined  earlier  and  other  works 

In  the  literature  is  that  only  one  real  random  variable  z(x)  Is  present  In 

both  the. -real  and  Imaginary  part  of  E(x,u)  whereas  In  the  calculation  of  P(I) 

13  14 

from  random  signals  * and  P(I)  for  light  scattered  from  a very  rough  sur- 
2-3 

face,  the  real  and  Imaginary  parts  of  the  amplitude  are  assumed  to  be  In- 
dependent Gaussian  random  variables.  Roughly  speaking,  P(I)  for  light  re- 
flected from  a moderately  rough  surface  Is  generated  by  real  Gaussian  variables 
z(x)  Instead  of  circular  complex  Gaussian  variables. 

z(x)  Is  expanded  In  terms  of  othogonal  functions  ^(x)  and  random  variables 

S1 

*(x)  - I Sjtj(x)  . (3) 

It  can  be  easily  verified  that  ^(x)  Is  a solution  of  the  Integral  equation 

J<z(x)z(*’)>  ^(x')  d2x  - X^fx)  (4) 

where  Xj  Is  the  eigenvalue.  Since  Sj  are  linear  combinations  of  Gaussian 
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random  variable  z^,  then 


< Sj  > ■ u (5) 

From  Eq.  (2)  and 

< > • Sij.X1  , (6) 

then  have  the  distribution  function  given  by  Eqs.  (5)  - (6),  or  , 

0(S1)dSi  * dSj  exp  (-  $2  / 2Xj)  /nJEttXj  . (7) 


The  Intensity  received  In  the  detector  I Is  explicitly  given 

2 2 

correct  to  the  second  order  In  k o by 

I..*  \ [(f0(7.5)  f0(*',D  ♦ 2 kV  f#(*.«>  f2(x'.u'j) 

♦ 2 f0(i.u)  f,(x'  .ff)  kz(x')  * k2  f^(x«u)z(x)  f,(x',u)  z(x')J 

d2x  d2x'  d2u  , (8) 

where  k2z2  Is  set  to  be  k2o2,  since  k2(z2-o2)  * 6 may  be  neglected. 
This  Is  justified  by  the  observation  that  the  average  of  S and  any 
function  W(kz) 

m 

< 6W(kz)  > ■ < «(kz)n>/n!  --  (9) 

" ' d(kz)n  kz-0 

Is  of  order  (kz)4,  since  n ■ 0 and  n ■ 1 In  Eq.  (9)  both  vanish. 
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Consequently  6 can  be  neglected.  Equation  (8)  can  be  written  as 
It  ■ I + I,  where  1^  Is  the  total  intensity  received  at  u and 
I is  from  the  Imaginary  part  of  the  amplitude  similar  to  Eq.  (8). 
Equation  (8)  can  be  written  In  the  form, 

1 - l0  + * I2  (10) 

where  Ij  contains  z(x)  of  1^  order.  I Is  defined  as  a sum  of 

terms  ? , I,,  and  L similar  to  Eq.  (10). 
o I i 

Equation  (3)  expresses  I as  a quadratic  function  of 

I ■ o ♦ E eiSi  + Z Y1jSiSj  (11) 

where  a , ^ , and  are  constants  obtained  after  Integration  of 

the  product  of  ^ and  various  functions  In  Eq.  (10). 

P(I)  Is  defined  by 

P(I)  • j ||  exp  (1  ItZ)  < exp  (-1  ItZ)  > (12) 

where  » I ♦ I,  and  < > Is  understood  to  be 

< ( ) > • nJdSt  DfS,)  ( ) . (13) 

The  calculation  of  P(I)  In  Eq.  (12)  can  be  surmarlzed  as  follows: 
first  solve  for  Eq.  (4)  to  obtain  ^ and  Sj.  Next  determine  a , 
and  Yjj  In  Eq.  (11).  Finally,  perform  the  Integral  P(I)  de- 
fined In  Eq.  (12).  Although  Eq.  (12)  looks  much  more  complicated  than 
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that  derived  for  a very  rough  surface,  the  calculation  in  general  can  be 
greatly  simplified  as  shown  later.  It  should  be  clear,  however,  that  for 
P(I)  for  a moderately  rough  surface,  the  real  and  Imaginary  of  the  amplitude 
cannot  be  separated  as  Independent  variables.  Because  of  this  reason, 

In  Eq.  (11)  are  not  diagonal. 


III.  SIMPLE  GEOMETRY 

In  order  to  demonstrate  how  P(I)  depends  on  various  parameters,  a 
simpler  scattering  geometry  for  numerical  calculation  Is  considered  here. 
It  should  be  clear  that  other  calculations  are  sufficiently  general  for 
many  similar  systems. 

We  consider  the  light  source  of  the  Gaussian  distribution. 


E0<*>  ' h 


exp  [-(x2  + y2)  /2T2]  Jlvl1 


reflected  from  a rough  surface  as  shown  In  Figure  1.  A lens  of  focal 
length  f Is  used  to  simplify  the  calculation.  E(x,u)  Is  explicitly  given 
by 

E(x,u)  • EQ(x)  G0(u)  xf  • exp  (l(lc-n).x)  • 

[l  + 1(kz  - n2)  . z - (k2  - n2)2  zZ/l\  (15) 

where  GQ(u)  « exp  (1  - (u2  + v2)j  ; and  k and  n are  the  propagation 

and  the  scattered  direction.  Notice  that  n > ku/f  depends  on  u . 

The  Intensity  I received  In  the  detector  Is 

I - Jd2u  Jd2xd2x'  Efl(x)  E0(x')  exp  1(k  - n)(x  - x') 

X [l  - (k2  - n2)2  o2  + 1(k2  - n2)(z(x)  - z(x'))+  (k2-n2)2z(x)z(x')J 
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or 


1 " Jd2u(j0(u.u)  + J2(u,u))  = I0  + I2  07) 

where  JQ  and  J2  are  functions  of  u only  and  J2  contains  the  random 
variables  z(x)z(x').  The  term  linear  In  z(x)  vanishes  as  shown  later. 

In  order  to  compute  P(I)  in  this  special  case,  we  find  It  simpler  to  follow 
many  previous  works  and  define  the  mutual  intensity  function  J2(u,u')  and 
solve  the  integral  equation 

Jo2(u,u')  ^(u‘)  d2u'  =*  A^tu)  (18) 

where 

J2(u,u')  * j^kz"nz^2  z(x)z(x')  exP  ) - 1(nx  - In'.x')] 

X E0(x)  E0(x')  d2xd2x'  (19) 

which  can  be  written  as 

j2(u,u’)  - t2  exp  £t2  (k  - -j-j-  j - T2  (tT-)  ] (20) 

where  u+  • u + u*  and  t < < T has  been  used.  Notice  that  J2 
becomes  a product  of  two  functions  of  u+  andu_,  respectively.  The  u_ 
dependence  comes  from  the  distribution  of  the  Incident  Intensity,  whereas 
u+  dependence  contains  the  surface  parameters.  This  Is  the  expected  re- 
sult for  any  Incident  distribution  much  coarser  than  t . In  the  limit 
of  small  t or  zero  correlation  length,  < J2  > depends  only  on  u_  , 
which  Is  the  Van-CIttel-Zenlke  theorem.^6 
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The  solution  from  Eq.  (18)  forms  a complete,  othonormal  set, 
from  which  we  expand 


$<Vnz>*<*>  •Kpfrt-")  * x]e(x)  d2x  -s  S1<*>1  (u) . (21) 

Similar  methods  used  earlier  show  that  the  random  variables  satisfy 


< Sj  > * 0 » 


<SlSj>  * *1  • 


(22) 

(23) 


Notice  that  the  expansion  Eq.  (21)  Is  different  from  Eq.  (3),  which 
Is  used  for  the  purpose  of  expressing  both  real  and  Imaglnery  part  of  the 
amplitude  In  terms  of  ip  . In  our  simple  geometry,  Eq.  (21)  Is  sufficient 
to  diagonalize  I2 


ij  -t  v 


(24) 


Although  In  principle,  we  still  have  to  find  z(x)  as  an  explicit  function 
In  order  to  express  the  linear  tern  In  Eq.  (17)  which  can  be  done  by  using 
Eq.  (21)  , In  our  practical  solution  the  linear  term  It- proportional  to  the 


Imaginary  part  of  l and  vanishes  since  and  are  real  as  a re- 
sult of  the  approximation  In  Eq.  (20). 


NMPIMi 


P(I)  from  Eq.  (12)  is 


-h 

9(1)  - j f=-  c°s  fz(I-I0)  - Ee^Jn  (1  + 4X2  Z2)  (25) 

•« 

where  * tan'1  (2x^z)  . 

We  will  use  the  one-dimensional  model  in  the  numerical  calculation 
to  discuss  how  P(I)  depends  on  various  parameters.  Eq.  (18)  then  Is  written 
In  the  form  for  the  normal  Incident  wave, 

1 

A o2k2c  j exp  (-Sc  u2  - c u2)  du*  <^(u)  * X^+^u)  (26) 


where  S * t/T  and  u is  the  unit  of  i , where  l Is  a dimension  of 

2 2 

the  detector,  c * ir(Tl/XT)  and  c/ir  may  be  Interpreted  as  the  number 
of  speckles  in  the  detector.  A is  a constant,  and  is  set  to  be  one,  since 
P(I)  depends  only  on  the  relative  magnitude  of  JQ  and  Jgi  and  IQ  has  been 
chosen  as  a parameter  In  the  following  discussion. 

The  effect  of  S < < 1 in  Eq.  (26)  In  general  Is  to  reduce  the  eigen- 
values of  x^  and  hence  reduce  the  contrast  as  easily  seen  later.  We  will 
discuss  first  the  case  S » 0 In  the  perturbation  calculation.  When  c 

■ 0,  there  Is  only  one  eigenfunction,  ^ * const.,  and  one  eigenvalue, 

2 2 2 2 2 
2 o k‘c.  For  small  c , the  kernel  can  be  expanded  on  2 o k c (1  • cu_) 

and  the  eigenfunctions  can  be  written  as  <^(u)  , 


tj(uJ  ■ Ai  ♦ 8jU  + C^u*  . (2 

By  substitution  of  Eq.  (27)  to  Eq.  (26),  we  Identify  the  terms  « un 
(n  ■ 0,1,2)  In  both  sides,  and  obtain  the  equation  for  eigenvalues  X1 
(1- 0.1,2) 
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0 


= o 


(28) 


- 2/3  c) 


(*i 


If  we  neglect  the  term,  and 


- 2/3  cj 
\Jl  c k2c2  . 


For  large  values  of  c the  kernel  becomes  a delta  function.  Many 
of  the  eigenvalues  are  approximately  equal  and  satisfy  the  same  rule 

Xj  » .2  c o2k2  . (29) 

In  Figure  2,  the  distribution  of  the  eigenvalues  Is  given.  It  Is  Interest 
Ing  to  note  that  there  Is  considerable  qualitative  differences  between  the 
distribution  of  the  eigenvalues  In  the  kernel  In  Eq.  (26)  and  the  sine 
function  used  In  Reference  5. 

Fpr  a set  of  equal  eigenvalues,  P(I)  can  be  easily  evaluated  by*6 

it  ] <«  «■>  [- 1 - v]"  ■ ijcrti"'  «*p 

r (n) 


where  o It  an  arbitrary  constant  with  R#n  > 0.  The  contrast  ct  on 

the  basis  of  Iq.  (90)  Is  given  by 


(30) 


IV.  NUMERICAL  RESULTS  AND  DISCUSSIONS 


P(I)  In  this  work  depends  on  three  parameters,  IQ  , c , and  S in 
Eqs.  (17)  and  (26).  The  dependence  on  IQ  only  shifts  the  origin  of 
P(I)  and  changes  the  value  of  c^  according  to  Eq.  (31).  The  depen- 
dence of  P(I)  on  c and  S Is  through  the  distribution  of  the  eigen- 
values shown  In  Figure  2,  which  then  effect  the  general  features  of  P(I) 
Let  us  examine  the  c dependence  by  setting  S « 0,  IQ  * 0.  In 
Figure  3,  we  plot  P(I)  for  c * 0.1,  1,  and  10.  c^  of  each  case  are, 
respectively,  1.24,  0.95,  and  0.53.  Notice  that  c^  can  be  greater  than 
one.  In  general  these  curves  behave  similarly  to  those  for  P(I,G)  for  a 
very  rough  surface.5  Recall  that  P(I,G)  Is  obtained  by  solving  for 
similar  eigenvalue  equations,  Eq.  (18),  except  the  power  in  the  integral 
In  Eq.  (25)  Is  replaced  by  one.  Instead  of  1/2.  It  Is  clear  from  Figure 
3,  th*t  for  small  c , P(I)  and  P(I,G)  are  very  different  at  small 
value  of  c ; however,  P( I)  and  P(I,G)  are  very  similar,  as  shown  in 
Figure  4,  where  c » 20  Is  used.  As  a matter  of  fact,  P(I)  behaves  very 

much  like  P(I,G)  at  a smaller  value  of  c . Eq.  (20)  has  a set  of  eigen 

values  of  about  the  same  order  of  magnitude.  We  also  have  fit  these 

curves  by  using  Eq.  (23),  where  IQ  and  n are  determined  from  the 

actual  curve  n * 2.7  and  6.5,  respectively,  for  P(I)  and  P(I,G).  The 
fact  that  Eq.  (23)  fits  P(I,G)  almost  perfectly,  has  also  been  demon- 
strated In  Reference  5,  where  the  fit  for  P(I,G)  of  smaller  value  of 
c was  also  shown  numerically  not  very  good. 


Next,  we  discuss  how  the  P(I)  depends  on  S , which  is  essentially 
the  reciprocal  of  the  number  of  correction  cells  in  the  target.  We 
should  emphasize  that  the  S dependence  is  very  much  related  to  the 
surface  model  and  the  scattering  geometry.  The  appearance  of  S in 
Eq.  (20)  also  explains  the  effect  of  S on  the  eigenvalues  in  Figure  2. 
First,  S does  not  change  the  eigenvalue  for  small  c . As  c becomes 
larger,  in  Eq.  (22)  is  reduced  according  to  the  scale  in  1//T".  As 
a result,  the  distribution  of  for  finite  S is  similar  to  that  of 

Sj  for  zero  S value  and  smaller  c . This  feature  of  the  effect  of 
finite  value  of  S (finite  number  of  the  correlation  cells)  is  shown  in 
Figure  5,  where  P(I)  is  plotted  for  c = 40  and  S = 0,  0.1,  and  0.2;  c^.  = 
.44,  .77,  and  .97,  respectively.  A comparison  of  Figures  5 and  3 shows 
that  P( I)  (C  * 40,  S * 0.2)  Is  rather  close  to  P (I ) (for  c * 1,  S = 0). 
Notice  that  for  S * 1,  there  will  be  only  one  eigenvalue  and  y (u)  - 
exp  (-  2 cy  ).  The  result,  however.  Is  not  meaningful  as  T > > t , 

S < < 0 has  already  been  assumed. 

The  only  detailed  experiment  relevant  to  our  calculation  is  in  Ref- 
erence 10,  where  P(I)  is  only  measured  for  c z 0.  P(I)  of  a point  de- 
tector Is  more  appropriate  for  direct  comparison,  although  the  qualitative 
dependence  of  P(I)  on  IQ,  and  S In  Reference  10  Is  In  agreement  with  the 
calculation.  It  Is  Important  to  note  that  IQ  may  be  a complicated  func- 
tion of  S , and  c10  , and  IQ  may  dominate  the  magnitude  of  ct  In  Eq. 
<24).  Thus,  the  experimental  control  of  IQ  Is  crucial  for  a qualitative 
comparison  with  the  calculation. 
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In  summary,  we  have  obtained  P(I)  for  a reflected  light  from  a 
moderately  rough  surface  as  a function  of  the  number  of  correlation 
cells  In  the  object,  and  the  size  of  the  scanning  aperture.  The 
calculation  supplements  Seribot's  work5  on  P(I)  for  a rough  surface. 


fl 
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FIGURE  CAPTIONS 


Figure  1. 
Figure  2. 

Figure  3. 

Figure  4. 

Figure  5. 


The  geometry  of  the  object  and  observation  plane. 

Distribution  of  the  eigenvalues  for  various  parameters.  Curves 
A,  B,  and  C are  for  S = 0,  c = 1,  10,  and  100,  respectively. 

B'  and  C'  are  for  S * 0.1,  c = 10  and  100,  respectively. 

P(I)  for  various  c and  S = 0.  Curves  A,  B,  and  C are  for 
c - 10,  1,  and  0.1,  respectively.  Each  curve  is  multiplied 
by  a factor  Y , YA  = 0.2,  Yg  = 1.3,  Yc  * 5.3 
P(I)  for  different  statistics.  Solid  curves  A and  B are 
based  non-circular  and  circular  Gaussian  statistics.  The 
dashed  curve  is  calculated  from  Eq.  (23). 

P(I)  for  c 3 40  and  various  S . Curves  A,  B,  and  C are  for 
S * 0,  0.1,  and  0.2,  respectively.  P(I)  is  multiplied  by  40. 
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